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Abstract 



o 



'Vibrio harveyi and Vibrio cholerae have quorum sensing pathways with similar design and highly homologous components in- 
cluding multiple small RNAs (sRNAs). However, the associated luminescence phenotypes of strains with sRNA deletions differ 
dramatically: in V. harveyi, the sRNAs act additively; however, in V. cholerae, the sRNAs act redundantly. Furthermore, there are 
striking differences in the luminescence phenotypes for different pathway mutants in V. harveyi and V. cholerae. However these 
differences have not been connected with the observed differences for the sRNA deletion strains in these bacteria. In this work, we 
present a model for quorum sensing induced luminescence phenotypes focusing on the interactions of multiple sRNAs with target 
mRNA. Within our model, we find that one key parameter - the fold-change in protein concentration necessary for luminescence 
Activation - can control whether the sRNAs appear to act additively or redundantly. For specific parameter choices, we find that 
differences in this key parameter can also explain hitherto unconnected luminescence phenotypes differences for various pathway 
^— ( mutants in V. harveyi and V. cholerae. The model can thus provide a unifying explanation for observed differences in luminescence 
^> phenotypes and can also be used to make testable predictions for future experiments. 
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\ Bacterial survival is critically dependent on regulatory net- 
works which monitor and respond to environmental fluctua- 
tions. An important example of such a regulatory network is 
the pathway responsible for bacterial "quorum sensing", com- 
monly defined as the regulation of gene expression in response 
to cell density |6]. Quorum sensing bacteria produce, se- 
crete and detect signalling molecules called autoinducers (AIs) 
which accumulate in the surroundings as the cell population in- 
creases. Differential expression of certain sets of genes occurs 
.when the local concentration of AIs exceeds a critical thresh- 
old. Several processes critical to bacterial colonization and vir- 
ulence e.g. biofilm formation, bioluminescence, and secretion 
of virulence factors 13211 were shown to be reg- 

ulated in this manner, leading to increased interest in character- 
izing quorum sensing based regulation in bacteria. 

The quorum sensing networks in Vibrio harveyi and Vibrio 
cholerae were recently analyzed in considerable detail ll26ll . 
The basic network components are highly homologous in the 
two species, to the extent that the bioluminescence genes from 
V. harveyi were used to characterize the re gula tory network 



in V. cholera using luminescence assays 11131 Il4 121[1 . In both 



species, the central regulatory module consists of multiple quo- 
rum regulatory small RNAs (qrrl-4 in V. cholerae and qrrl-5 



'Corresponding author 

Email address: kulkarniOvt . edu (Rahul V Kulkami) 
' The corresponding bioluminescence genes are absent in V. cholerae. 



in V. harveyi) which control levels of the master regulator for 
quorum sensing: LuxR in V. harveyi and HapR in V. cholerae. 
LuxR/HapR levels are maintained below the threshold for lu- 
minescence activation at low cell densities due to repression by 
the small RNAs (sRNAs), whereas at high cell densities quorum 
sensing leads to a reduction in sRNA production rates, thereby 
increasing LuxR/HapR levels above the threshold leading to lu- 
minescence activation. By observing luminescence levels as 
a function of cell density for different mutants (corresponding 
to different deletions in the quorum sensing pathway compo- 
nents), several characteristics of pathway structure and function 
were inferred. 

The above studies documented striking differences in lumi- 
nescence phenotypes in the two species even though the regu- 
latory components of the pathways are very similar. The most 
dramatic differences were seen in the luminescence phenotypes 
of the qrr sRNA mutants. In V. cholerae, the four qrr sR- 
NAs acted redundantly 1,141 - all mutants with only one sRNA 
present had luminescence phenotypes that were identical to the 
WT luminescence phenotype. In contrast, the corresponding 
sRNAs in V. harveyi acted additively such that different mu- 
tants with only one sRNA present had distinct luminescence 
phenotypes compared to the WT phenotype. Thus in V. harveyi 
all the sRNAs must be present in order to mimic the wild-type 
luminescence phenotype [30]. Apart from these differences in 
the luminescence phenotypes of the sRNA mutants, there were 
also significant differences in the luminescence phenotypes of 
strains corresponding to deletions of upstream pathway ele- 
ments in the two species. An important challenge for compu- 
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tational analysis of quorum sensing pathways is to present a 
unifying explanation for the various, apparently unrelated, dif- 
ferences in the luminescence phenotypes for the two species 
despite the fact that the pathway elements are very similar (see 
figure [U. 
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Figure 1: A schematic of the quorum sensing gene networks. (Left) V. harveyi 
and (Right) V. cholerae employ multiple AIs whose signals are integrated to- 
gether in order to regulate either LuxR or HapR. V. harx'eyi produces and moni- 
tors the concentrations of three dilferent AIs (HAI- 1 , CAI- 1 , and AI-2), while V. 
cholerae produces and monitors the concentrations of two different AIs (CAI-1 
and AI-2). Via very similar' phosphorelay networks composed of highly homol- 
ogous components, the bacteria transduce the signal produced by the external 
AI concentrations through the network. In both bacteria, the sensors transfer 
phosphate groups to the protein LuxU when the external concentration of AIs 
is low. LuxU then passes the phosphate groups to the protein LuxO which, 
when phosphorylated, is responsible for the production of sRNAs. The flow of 
phosphate groups slows and then reverses when the external concentration of 
AIs continues to increase, thus reducing the production of the sRNAs. 

Previous work on modeling quorum sensing globally has de- 
veloped a framework for measuring a bacterium's ability to 
sense its microenvironment [27] ■ By modeling the relative con- 
centration of autoinducers inside and outside the cell, quorum 
sensing bacteria can be characterized based on a 'sensing poten- 
tial' and its relation to the associated activation properties, e.g. 
the critical population density ll27ll . Modeling has also been 
done at the genetic interaction level, in particular the interac- 
tion of sRNAs with target mRNAs. The solutions for the corre- 
sponding rate equations (for a range of parameter values) con- 
tained a sharp transition from a steady state wherein the target 
mRNA was strongly repressed to one in which the sRNA was 
strongly repressed IUImEIIIIII. Since the WT lumines- 
cence phenotype also showed a sharp transition as cell density 
increased, it was initially suggested that this transition corre- 
sponded to the sharp transition seen in the sRNA-target rate 
equations. However, recent experimental results provide indi- 
cations that this is not necessarily valid and correspondingly the 
picture needs to be revised. 

Experiments in V. cholerae showed that the expression lev- 
els of the virulence regulator AphA ll 1 111 to be about three-fold 
lower in WT at low cell densities compared to a AhapR mu- 
tant. This indicates that WT V. cholerae maintains HapR at low 



but significant levels at low cell densities (such that it can ef- 
fectively repress AphA to the extent noted) rather than fully re- 
pressing it. Furthermore, experiments in V. harveyi examining 
regulation of additional targets by LuxR indicated that LuxR 
levels change in a graded manner as opposed to a sharp, ultra- 
sensitive switch 11331]. Thus, there is a need for computational 
analysis of sRNA-target regulatory interactions in the context 
of the quorum sensing pathway, which is consistent with these 
experimental results and which also provides a unifying expla- 
nation for observed luminescence phenotypes. 

In what follows, we will present a simplified model for lu- 
minescence regulation during quorum sensing in V. harveyi and 
V. cholerae, which is an extension of work previously done by 
the authors of this paper |2]. For a given choice of parameters, 
the model accounts for the dramatic differences in the lumi- 
nescence phenotypes for the sRNA mutants in the two species 
based on a single parameter difference. The analysis also pro- 
vides a unifying explanation for currently unrelated differences 
between the luminescence phenotypes of different mutants in 
the quorum sensing pathways and gives rise to testable predic- 
tions for future experiments. This work thus provides a frame- 
work for systems-level analysis of the quorum sensing pathway 
in the V. harveyi and V. cholerae while complementing previous 
models of V. fischeri 11121 124, 12511 and suggests future experi- 
ments that can help in further unravehng the function of this 
critical regulatory pathway. 

1.1. Overview of experimental results 

We begin with an overview of the two pathways and associ- 
ated luminescence phenotypes in the two species. A schematic 
representation of the two pathways is shown in figure [T] The 
core elements are the same in both species: a multi-component 
phosphorelay involving sensor proteins (which can function as 
kinases as well as phosphatases), the phosphotransfer protein 
LuxU, and the response regulator protein LuxO. Phosphory- 
lated LuxO is responsible for the activation of multiple qrr sR- 
NAs which in turn repress the quorum sensing master regulator 
(LuxR in V. harveyi and HapR in V. cholerae). 

The pathways do exhibit some differences in the number of 
autoinducer synthase/sensor protein pairs and in the number of 
sRNAs present. V. harveyi has three known autoinducer syn- 
thase/sensor protein pairs whereas V. cholerae has only two 
known autoinducer synthase/sensor protein pairs. Furthermore, 
V. harveyi has five qrr sRNAs as opposed to four in V. cholerae 
fl^. However, our current understanding indicates that these 
differences are not significant under the conditions tested. For 
example, it was shown that qrrS in V. harveyi is not quorum 
sensing regulated or expressed under normal conditions 130(1 
and one of the autoinducer synthase/sensor protein pairs in V. 
harveyi has minimal effects on quorum sensing based regula- 
tion [9*]. Thus, both pathways can effectively be considered 
as having two autoinducer synthase/sensor protein pairs and 
four qrr sRNAs. Furthermore, the pathway components are 
highly homologous, e.g. LuxR is greater than 90% identical 
to HapR. However, despite these common features and similar- 
ities between components, the luminescence phenotypes show 
dramatic differences as detailed below. 
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Figure 2: Experimental Wild-Type luminescence curves for V. harveyi and V. 
cholerae. The solid, red curve represents the change in luminescence relative 
to optical density (OD) for V. han'eyi. There is a smooth transition in lumines- 
cence near OD 10"' as the distribution of cells switch from "off' to "on" t30ll . 
The dashed, red curve represents the change in luminescence relative to OD for 
V. cholerae. There is a sharp transition in luminescence near OD 10° as the dis- 
tribution of cells switch from "off' to "on" 1 14]. The vertical solid and dashed 
lines represent possible OD concentrations that correspond to the beginning of 
the cells in the population reaching a LuxR/HapR concentration necessary for 
luminescence for V. han'eyi and V. cholerae respectively. Regions indicated by 
(1), (2), and (3) reflect the relative protein distributions labeled similarly and 
shown in figure |4] 



While the most dramatic differences in luminescence phe- 
notypes occur in the qrr and luxU mutants, distinct qualita- 
tive differences exist even in the WT phenotypes, see figure |2] 
These qualitative differences serve as motivation for our mod- 
eling and help in discerning a key difference between the two 
genetic networks that can potentially explain the various differ- 
ences in luminescence phenotypes. The luminescence curves 
for WT strains of V. harveyi and V. cholerae (based on experi- 
mental data from llsoll and |14]) are shown in figure |2l In both 
cases, the luminescence per cell begin at a high value since the 
initial state corresponded to a dilution of the high cell density 
culture which was maximally bright. As the colony density in- 
creases, the luminescence level drops until a critical cell density 
is reached, after this critical point there is a subsequent rise in 
luminescence back to the initial level. While the luminescence 
curves of WT V. harveyi and V. cholerae look similar, there 
are important differences between the two curves. Wild-type V. 
harveyi showed an almost symmetric parabola centered around 
ODgoo ~ 0.1 [30]; however, wild-type V. cholerae showed a 
continued decline in relative light unit (RLU) output until the 
colony reached an ODsoo ~ 10. The luminescence levels then 
increased by several orders of magnitude over a timescale dur- 
ing which cell density changed by a small factor A fold) 114 1 . 

The luminescence phenotypes of strains corresponding to 
deletions of various pathway elements also depicted important 
differences between the two species. As mentioned in the In- 
troduction, luminescence curves of qrr sRNA mutants in the 
two species suggested that the sRNAs functioned additively in 
V. harveyi fid\ but were redundant in V. cholerae fTT). Another 
striking difference was seen in the luxU mutant which was al- 
ways bright regardless of cell density in V. harveyi whereas the 
luxU mutant showed a density-dependent luminescence pheno- 



type in V. cholerae. Furthermore, while deletion of the sensor 
kinases (e.g. for the cqsSJuxQ mutant) changed the lumines- 
cence phenotype with respect to WT for V. harveyi, the corre- 
sponding WT and deletion mutant luminescence curves were 
almost identical for V. cholerae 11311 . These observations based 
on experimental luminescence curves lead to some important 
questions which need to be addressed: 

1) How can we understand changes in RLU (Relative Light 
Unit)/cell over several orders of magnitude corresponding to 
small changes in cell density? 2) How are the phenotypes dra- 
matically different despite the basic components/circuitry being 
the same? 3) Is there a unifying explanation for the seemingly 
unrelated differences in luminescence phenotypes for different 
mutant strains? 



2. Methods 

2.1. Modeling framework 

In order to address the issues raised above, we will first 
discuss the modeling framework and key assumptions of our 
model. They are schematically illustrated in this section and 
more quantitatively developed in following sections. 

We assume that the measured luminescence levels per cell 
are proportional to the rate of transcription of the luminescence 
genes. Since these genes are activated by the quorum sensing 
master regulators, their transcription rate is a function of cel- 
lular concentrations of LuxR/HapR. We assume that this func- 
tion has a sharp threshold; as a simplification we represent it 
by a step function such that cells with LuxR/HapR concentra- 
tions below the threshold produce no light whereas cells with 
LuxR/HapR concentrations above the threshold produce maxi- 
mal luminescence. 

Since the experimentally measured quantity is the popula- 
tion average of the luminescence output/ceU, we need to con- 
sider the steady state distribution of LuxR/HapR levels across 
all cells. Recent work showed that the steady state protein dis- 
tribution for proteins can be characterized as a Gamma distribu- 
tion |4]. Accordingly, we represent the LuxR/HapR distribution 
by a Gamma distribution with a given variance and whose mean 
value is determined by solving the rate equations of our model 
(see next section). 

With the assumptions mentioned above, we can make signif- 
icant inferences about quorum sensing networks based on the 
luminescence data. The change in RLU/cell over several or- 
ders of magnitude corresponds to the steady state distribution 
for LuxR/HapR crossing the luminescence activation thresh- 
old (see figure[3]l. Thus the mean concentration of LuxR/HapR 
must change by the minimal amount indicated in the figure dur- 
ing the transition from the 'dark' phenotype to the maximally 
luminescent phenotype. The WT luminescence curves indicate 
that this change occurs gradually in V. harveyi (positions (1) 
and (2) in figure |2]) as compared to V. cholerae (positions (2) 
and (3) in figure|2]i. Since the change in mean HapR levels in V. 
cholera (at OD ~ 1 .0) occurs without a corresponding signifi- 
cant change in cell density, it is unlikely to be driven solely by 
quorum sensing. Instead we infer, based on the luminescence 
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Figure 3: An illustration depicting luminescence activation as LuxR/HapR con- 
centrations cross a sharp threshold for activation. More (less) than the thresh- 
old, luminescence is (not) activated. The dilferent distributions depicted in this 
illustration are examples of LuxR/HapR concentration distributions and the cor- 
responding luminescence profile for a few examples, (left - cyan) A protein 
distribution that remains below the threshold regardless of cell density, (mid- 
dle - green) A protein distribution that transitions across the threshold and is a 
function of cell density, (right - magenta) A protein distribution that is entirely 
past the threshold regardless of cell density. 



phenotype, that there is a sharp rise in HapR levels around OD 
~ 1.0 in V.cholerae. One potential cause for this rise is a fur- 
ther reduction in the available regulatory sRNAs allowing for 
more available hapR transcripts. A possible source of the sRNA 
reduction is that as the cells move into stationary phase from 
growth phase, there is a decrease in the production of the Hfq 
chaperone [10]. A decrease in Hfq corresponds to a decrease in 
the concentration of sRNA-Hfq complexes which are necessary 
to regulate the target mRNA. Recent experiments in V. choleme 
have indeed found evidence for a sharp rise in HapR levels at 
OD~ 1.0 ['28']. 

In contrast, the transition in the WT luminescence phenotype 
for V. harveyi occurs at lower OD values and is more gradual 
suggesting that it is driven by the quorum sensing pathway. This 
observation leads to the suggestion that the crucial difference 
between the two species Ues in the location of the threshold for 
luminescence activation: in V. harveyi, quorum sensing based 
regulation suffices for moving the steady state LuxR distribu- 
tion across the threshold, whereas in V. cholerae this requires 
an additional jump in HapR levels at OD ~ 1.0. 

2.2. A minimal model for luminescence activation 

We focus on quorum sensing pathway elements correspond- 
ing to the production of sRNAs, the transcription of the target 
mRNA (luxR or hapR), and the interaction between the sRNAs 
and target mRNA. We start with a model containing only one 
sRNA species and neglect autoregulation of the target protein. 
Then we add the contributions of multiple sRNAs and autoreg- 
ulation to the model. 

The basic equations for a simplified model of sRNA-target 
interaction have been introduced and analyzed in previous work 
il4- 16, 23] and are given below (equations ([T]i and (|2|i). Con- 
sider first the case of a single sRNA species regulating one 
target mRNA species. If [x] denotes the concentration of the 



sRNA and [y] the concentration of the target mRNA, the corre 
sponding equations are: 
d[x] 



—J- = - [y] - px [x\ , 
dt 

^] 

dt 



—— = ky - y{x\ [y] - py [y] , 



(1) 



(2) 



where the ^'s are the production rates of each species, the p's 
are the degradation rates of each species, and y is an effective 
parameter for mutual degradation of sRNA and target mRNA. 

To generalize the above equations ([TJ and (|2]i while taking 
care of the effective parameter constraints (see Appendix A), we 
include the effects of 1 ) multiple sRNAs regulating luxR/hapR 
and 2) autoregulation of LuxR/HapR Ol IlTtl . The correspond- 
ing equations are, 
d[xi] 



dt 
d[y] 
dt 



= kx, - 7i[xi][y] - Px,[xi], 

JCv 



(3) 



i + ily/lyo]) 



J]i7i[xi]ly])-f^ybl (4) 



The constant [yo], represents the threshold concentration for 
binding of the target protein to its own mRNA. When the target 
protein is bound to the promoter region, transcription of the 
target gene is effectively blocked. 

Bioinformatic analysis il4ll indicates that the 32 bp region 
in the qrr sRNAs which is involved in regulation of hapR/luxR 
is absolutely conserved for all the sRNAs. Thus, we make the 
assumption that all the sRNAs have the same affinity for the 
target mRNA, i.e. we set ji - y. We further assume that the 
degradation rates of all sRNAs are the same {pi — p). However, 
the model does consider differences in the sRNAs production 
rates (A:,,) as demonstrated by experiment ifsT]. 

At steady state, the mean protein concentration is the mean 
mRNA concentration scaled by a constant - the ratio of the pro- 
tein translation rate to the protein degradation rate. Therefore, 
we use the scaled mRNA concentration in place of the protein 
concentration (see Appendix). 

To make the connection to luminescence curves, we have to 
consider the distribution of protein levels across cell popula- 
tions. Recent work by Friedman et al. showed the distribution 
of the protein concentration per cell for the colony can be repre- 
sented by a Gamma distribution [4]. Furthermore, recent flow 
cytometry work showed the distributions of fluorescence per 
cell from a luxR-gfp fusion had a nearly constant variance for 
a variety of conditions related to the concentration of AIs l33|] . 
Similarly, luminescence output from Vibrio harveyi was shown 
to be heterogenous across cell populations [1]. Therefore, we 
model the protein distribution as a Gamma distribution with a 
fixed variance. The mean of the distribution is obtained from 
the equations above for a given choice of parameters. Using 
this framework, we show in the following section how a single 
parameter difference can account for the vastly different lumi- 
nescence phenotypes of V. harveyi and of V. cholerae. 

3. Results and Discussion 

In this section we show how the minimal model discussed 
above with only one essential difference (the threshold for lu- 
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minescence activation) between the V. harveyi and V. choleme 
pathways can explain the observed differences in luminescence 
phenotypes as well as lead to testable predictions. 

We note that bacterial colonies are observed to change their 
luminescence production by many orders of magnitude in a rel- 
atively short amount of time, see figure|2] However, the changes 
in the level of the master regulator proteins and sRNAs are not 
nearly as dramatic ifsol 331. We interpret this as indicating that 



a significant fraction of all the cells in the colony reach the 
conditions necessary for luminescence activation upon a small 
change in the master regulator protein levels. We model this as 
corresponding to a significant fraction of the master regulator 
distribution moving across sharp threshold values of concentra- 
tions necessary to activate luminescence (see figure[3l). 

3.0.1. The protein distributions for WT strains 

Using the model equations with parameter values guided by 
experiment (see Appendix), we plot the the distribution of the 
protein concentration for a WT colony (representing either V. 
harveyi or V. cholerae) at the low-cell density limit, high-cell 
density limit, and entering stationary phase labeled as positions 
(1), (2), and (3) respectively, see figureH] Since the protein dis- 
tributions for WT and all the mutants in either V. harveyi or V. 
cholerae are not available, we plot the distributions with respect 
to fold changes relative to the mean protein concentration for a 
WT colony at the low-cell density limit0 Specifically, the first 
two distributions in figure |4] are representative of the maximal 
relative change in protein concentration in going from low-cell 
density to high-cell density based on changes due to quorum 
sensing alone. The third distribution in figure|4]is the resulting 
distribution after the final reduction in sRNA production lead- 
ing to a rise in HapR due to entering stationary phase. 

The distributions at positions (1) and (2) in figure |4] repre- 
sent the maximally dark and maximally bright WT V. harveyi 
colonies, respectively. Similarly, the distributions at positions 
(1-2) and (3) in figure |4]represent the maximally dark and max- 
imally bright WT V. cholerae colonies, respectively. Ideally 
none of the bacteria in the dark colony should be "on" and 
none of the bacteria in the bright colony should be "off', there- 
fore we set the threshold of light activation for V. harveyi at a 
fold change directly in between the two distributions - depicted 
as the solid, vertical line in figure H] Since experiments have 
shown that the activation of V. cholerae to occur at a larger cell 
density than V. harveyi, we propose the threshold of light activa- 
tion for V. cholerae to be at a larger fold change - depicted as the 
dashed, vertical line in figure|4] As indicated in the figure, this 
corresponds to luminescence activation occurring in V. harveyi 
using quorum sensing alone, whereas for V. cholerae lumines- 
cence activation requires both transition to the high-cell density 
limit for the quorum sensing pathway and additional changes 
in HapR levels associated with entry into stationary phase. In 
what follows, we will discuss how assuming V. cholerae has a 
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-At the time of this work, the experimental protocol for measuring the ex- 
act protein concentration per cell in vivo was not available and has now only 
recently been published |[29|]. 



Figure 4: The distributions of the protein concentration across a WT bacterial 
colony from the model for the: (1) low-cell density limit, (2) high-cell density 
limit, and (3) entering stationary-phase limit. The x-axis depicts the fold change 
difference relative to the mean protein concentration value for a WT colony 
at the low-cell density limit. The solid, vertical bar between distributions (1) 
and (2) and the dashed, vertical bar vertical between distributions (2) and (3) 
represent the threshold values for luminescence for V. harveyi and V. choleme, 
respectively. 



different threshold of light activation than V. harveyi can con- 
sistently explain the differences in the sRNAs and luxU mutant 
phenotypes. 

3.0.2. Additivity vs redundancy 

We account for each of the four active qrr sRNAs having 
a different production rate and set the rates with the follow- 
ing hierarchy: qrr4 > qrr2 > qrr3 > qrrl, which is consistent 
with experimental results in V. harveyi [3Q]. Figure|5]shows the 
distributions of the protein concentrations for mutant colonies 
containing only one of the four active qrr sRNAs for both V. 
harveyi and V. cholerae - each sRNA mutant is represented as 
a different shade of green in figure |5] 

In the low-cell density limit, position (1) in figure|5] the dis- 
tributions all have regions extending past the threshold for lu- 
minescence in V. harveyi. This is a representation of the qrr 
"additivity" response seen in V. harveyi as all qrrs are needed to 
prevent any appreciable region of the protein distribution from 
extending past the threshold in the low-cell density limit ll30ll . 
For the sRNA mutants, the regions of the distributions in the 
low-cell density limit that extend past the threshold represent 
the amount of bacteria in the colony that are "on" regardless of 
cell density. 

The story is difl'erent from the perspective of V. cholerae. 
In the high-cell density limit, position (2) in figure |5] the dis- 
tributions are all below the threshold for luminescence in V. 
cholerae, which corresponds to complete light repression and 
mimics the WT V. cholerae response iflitl . We suggest that 
once the final reduction in sRNA production occurs, e.g. en- 
tering stationary phase, all the distributions cross the threshold 
for luminescence, position (3) in figure|5] The resulting pheno- 
type looks to be the same as the WT V. cholerae phenotype with 
the conclusion that the sRNAs act "redundantly". However, the 
prediction from our model is that the sRNAs behave the same 
in both V. harveyi and V. cholerae, but the associated thresholds 
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for luminescence are different. 
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Figure 5: The distributions of the protein concentration across a mutant colony 
containing only one active qrr sRNA from our model at the: (1) low-cell density 
and (2) high-cell density limits. The x-axis depicts the fold change difference 
relative to the mean protein concentration value for a WT colony at the low- 
cell density limit. The solid, vertical bar between distributions (1) and (2) and 
the dashed, vertical bar vertical between distributions (2) and (3) represent the 
threshold values for luminescence for V. haneyi and V. cholerae respectively. 



3.0.3. The luxU mutant 

The luxU mutant is another example of a difference in lumi- 
nescence phenotypes between V. harveyi and V. cholerae. The 
protein LuxU is responsible for coupling the autoinducer input 
signal to the rest of the quorum sensing network, see figure[T] If 
LuxU is removed from the pathway, the total sRNA transcrip- 
tion rate would drop to minimal levels, and the system would 
no longer respond to changes in cell density. Therefore if the 
quorum sensing pathway is the only factor controlling the lumi- 
nescence phenotypes, removal of the the luxU gene should re- 
sult in a bright, density independent phenotype. For V. harveyi, 
this is indeed the case - the luxU mutant is bright regardless of 
cell density. 

The story, as before with the sRNAs, is different with V. 
cholerae. In V. cholerae, the luxU mutant shows a density de- 
pendent luminescence phenotype, but the shape of the lumi- 
nescence curve is different from the canonical quorum sensing 
luminescence curves 112 ill . In the low-cell density limit, there 
is a detectable level of light production that is larger than WT 
value or any of the sRNA mutants values but much less than the 
maximal level of light production. This low level of lumines- 
cence remains stable for a significant portion of the exponential 
phase, and then sharply increases to the maximum level of lumi- 
nescence - a feature present in most V. cholerae luminescence 
curves. 

Our model reproduces this observed luxU mutant behavior 
in V. harveyi and V. cholerae. In figure |6] there are only two 
distributions: one for the high-cell density limit (position (2)) 
and one for the high-cell density limit entering stationary phase 
(position (3)). From the perspective of LuxR/HapR regulation, 
the removal of luxU effectively decouples the quorum sens- 
ing pathway from the outside inputs. Therefore, the system 
effectively starts at the high-cell density limit, and the associ- 
ated protein distribution is always past the V. harveyi lumines- 
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Figure 6: The distributions of the protein concentration across a mutant colony 
where luxU has been removed from the system from our model at the: (2) high- 
cell density limits and (3) high-cell density limit entering stationary phase. The 
X-axis depicts the fold change difference relative to the mean protein concentra- 
tion value for a WT colony at the low-cell density limit. The solid, vertical bar 
between distributions ( 1 ) and (2) and the dashed, vertical bar vertical between 
distributions (2) and (3) represent the threshold values for luminescence for V. 
harveyi and V. cholerae respectively. 



cence threshold. This results in a fully bright, density indepen- 
dent phenotype, see figure |6] However, the distribution asso- 
ciated with the high-cell density is only partially across the V. 
cholerae luminescence threshold resulting in a small concentra- 
tion of the cells being "on" and a majority being "off'. The V. 
cholerae colony will remain in this state until it enters station- 
ary phase where the protein distribution completely crosses the 
V. cholerae luminescence threshold, see figure |6] 

Finally, the luminescent behavior of the cqsS and luxQ dou- 
ble mutant in V. cholerae is also consistent with the model. Es- 
sentially, this double mutant shows a WT response even though 
the both autoinducer sensors are removed \\y\. In our model, 
this would correspond to the system starting in the WT high- 
cell density limit, position (2) in figure |4] which is below the 
threshold for luminescence in V. cholerae. Therefore, the ob- 
served phenotype should be nearly identical to WT. 

By just having two thresholds separating three distinct re- 
gions of protein regulation, our model is able to consistently 
link the sRNAs acting additively in V. harveyi llsoll . the sRNAs 
acting redundantly in V. cholerae lll4ll . and the density depen- 
dent phenotype in V. cholerae for the luxU mutant [21]. With 
the relative positions of the thresholds and protein distributions 
now in place, we now discuss the predictions that come from 
our model. 



3.1. Predictions 

Current experimental techniques can produce a variety of dif- 
ferent mutant strains of V. harveyi and of V. cholerae. Depend- 
ing on the genes and sRNAs being removed from the strain, the 
experimental techniques generate even up to triple knock-out 
mutants (and possibly more if required). Since our simplified 
model, with the given choice of parameters, reproduces features 
of the observed luminescence phenotypes, it is of interest to ex- 
amine the model predictions for luminescence phenotypes of 
different gene and sRNA mutant strains that should be experi- 
mentally feasible to test. 
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The model distinguishes the varying behaviors of V. harveyi 
and of V. choleme as a difference in the threshold protein con- 
centration of the master regulatory gene, and the concentration 
of the master regulatory gene at any position is determined by 
the associated production rate of the sRNAs. Therefore, there 
is an effective total sRNA production rate that coincides with 
the distribution of the master regulatory protein being centered 
at a given threshold value. We refer to this critical value of total 
sRNA production as kc- 

Since we assume the threshold values are different for V. 
harveyi and V. choleme, their associated critical value of total 
sRNA production, kc, is different. Each mutant has an associ- 
ated total sRNA production rate at low-cell density and high- 
cell density limits. A hierarchy of sRNA production rates for 
different mutants and colony cell densities relative to kc ex- 
plains currently seen phenotypes, and we will use this hierarchy 
as a basis for predicting new phenotypes. 

WTi > sRNAi >kc> WTh > sRNAh > AU > AO, (5) 
WT, > sRNAi > WTh > sRNAh > AU > kc > AO. (6) 

Equations ^ and (|6]l represent the hierarchies for V. harveyi 
and V. cholerae, respectively. Those rates greater than kc cor- 
respond to "dark" phenotypes, and those rates less than kc cor- 
respond to "bright" phenotypes. In equations (|5]l and (|6]l, WTi 
and WTh represent the total sRNA production rate for wild-type 
bacteria in the low-cell density and high-cell density limits be- 
fore the transition to stationary phase. sRNAi and sRNA/, are 
the sRNA production rates for any mutant with at least one ac- 
tive sRNA removed from the system in the low-cell density and 
high-cell density limits before the transition to stationary phase. 
Finally, AU and AO in equations (|5]i and (|6|l are the sRNA pro- 
duction rates for the mutants where LuxU and LuxO has been 
deleted, respectively. Now that the hierarchy is established, we 
discuss below the resulting predictions. 

One way to explore the different quorum sensing responses 
of the network is to add an external concentration of autoinduc- 
ers to a low-cell density colony, also know as "cross-feeding". 
The additional autoinducers will "trick" a colony into behaving 
as if it is in the high-cell density limit which causes a transi- 
tion in the total sRNA production rate. Since the production 
rates specifically dependent on cell density, WTi, sRNAi, WTh, 
and sRNAh are separated by kc in equation (|5]l, the model pre- 
dicts a low-cell density colony of wild-type or any sRNA mu- 
tant V. harveyi will start to luminesce when extra autoinducers 
are added to the colony. 

However, for V. cholerae, the production rates specifically 
dependent on cell density are all greater than kc in equation (|6]l. 
Even the production rate of the luxU mutant is greater than kc- 
Therefore, the model predicts that a low-cell density colony of 
wild-type or any sRNA mutant V. cholerae will remain dark 
when extra autoinducers are added to the colony. Also, the 
model predicts this outcome for any mutant V. cholerae cor- 
responding to a total sRNA production rate greater than kc, in- 
cluding the luxU mutant and the cqsS and luxQ double mutant. 

The model also predicts cases where the sRNAs act "addi- 
tively" in the luxU mutant V. cholerae. Since the total sRNA 



production rate associated with the luxU mutant ^ is adjacent 
to kc, reducing the total sRNA production rate to a value less 
than kc will result in light production. This could be achieved 
by combining luxU with sRNA mutants. Therefore, the differ- 
ent sRNA triple mutants in combination with the luxU mutant 
for V. cholerae should show the associated HapR concentration 
changing in a graded manner Thus our model predicts that, in 
a luxU mutant background, the different triple sRNA mutants 
will appear to behave additively with regards to the lumines- 
cence - a phenotypic response similar to sRNA mutants in WT 
V. harveyi. 

3.2. Discussion 

In this study, we have shown how a simple set of equations, 
with appropriate choice of parameters, can effectively mimic 
the quorum sensing luminescence phenotypes of V. harveyi and 
V. cholerae. While the components of the quorum sensing reg- 
ulatory network in each of the bacteria are biologically similar 
in both homology and function, there are striking differences 
in luminescence phenotypes for the same mutant, e.g. luxU. 
Even the sRNAs, which are virtually identical in their sequence 
specificity to the target gene, act additively versus redundantly 
for V. harveyi and V. cholerae. 

We account for the striking differences by suggesting that 
the threshold concentration of the master protein needed for the 
bacteria to start luminesce activation is larger in V. cholerae than 
V. harveyi. The larger threshold concentration correspondingly 
implies the need for a mechanism that increases the levels of 
the master regulator in addition to the increases due to quorum 
sensing. The increase in master regulator levels can be effec- 
tively modeled as a sharp drop in sRNA productions rates and 
one possible source of this reduction can arise from the transi- 
tion from exponential growth phase to stationary phase. 

We considered solutions of the model equations for specific 
parameter choices motivated by experiments and analyzed the 
effect different mutants have on the sRNAs' production rates. 
In V. harveyi, the removal of either LuxO or LuxU causes a suf- 
ficient reduction in the sRNAs' production rate to result in the 
bacterial colony achieving maximal luminescence at any cell 
densities. Only the removal of LuxO from V. cholerae results 
in a similar response. Removing LuxU does not drop the sR- 
NAs' production rates enough for the bacteria to luminesce at 
any cell density. The extra reduction in the sRNAs' production 
rates from the transition to stationary-phase is required for the 
luxU mutant of V. cholerae to luminesce. 

We note that differences in the rate parameters between the 
two species could also account for some of the observed differ- 
ences in phenotypes. However, quantifying the multitude of re- 
action rates in both organisms is a challenging task experimen- 
tally. Also, even if significant differences in the rates between 
the two species are found, it is not clear if they will account 
for the observed dramatic differences in phenotypes. Instead 
of considering a multitude of parameter differences as the ex- 
planation for the observed phenotypes, our model suggests that 
changes in a single parameter can provide a unifying explana- 
tion for all the observed differences. Furthermore our model 
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leads to testable predictions which can easily be validated ex- 
perimentally (e.g. changes in luminescence upon crossfeeding). 

We further note that an alternative explanation for the dif- 
ferences in the luminescence phenotypes (in particular for the 
luxU mutant phenotypes) is based on the observation that the 
VarA/S-CsrA pathway interacts with the the quorum sensing 
pathway in V. cholerae \ 1 3i1 . While the corresponding interac- 
tion has not been studied in V. harveyi, it is likely that VarA/S- 
CsrA pathway interacts with the quorum sensing pathway in a 
similar fashion in V. harveyi. If, however, it turns out that the 
interaction between the two pathways is absent in V. harveyi, 
then this observation could account for some of the differences 
seen. In this case, the fold-change required for luminescence 
activation would be higher in V. cholerae due to repression by 
CsrA taking place in V. cholerae but not in V. harveyi. However, 
even with the deletion of csrA in V. cholerae, the luxU mutant 



sionless parameters; x - (}ixlkx)[x\, y - {jJiy I ky)[y], a = 
{yky)/{fi,xfiy), and yS - (y^.v)/(jUA/Uy) ™ equations (O and (|2]i so 
that the corresponding equations at steady state become: 



still shows a density-dependent phenotype in V. cholerae 11 1311 in 
contrast to the observed phenotype in V. harveyi. Thus there is a 
clear difference between the luxU mutant phenotypes between 
the two species which cannot be ascribed to potential differ- 
ences in the interactions between the quorum sensing pathways 
and the VarA/S-CsrA pathway. 

Within our model, the relationship between the threshold 
concentration and the total of all the sRNAs' production rates 
leads to experimental predictions. The first prediction is the 
inabiUty to prematurely initiate luminescence in a low-cell den- 
sity colony of V. cholerae through the addition of a large con- 
centration of autoinducers. Thus cross-feeding based activation 
of luminescence should work in V. harveyi but not in V. cholerae 
We also predict that a luxU mutant of V. cholerae combined 
with sRNAs mutants will result in a phenotype where the sR- 
NAs act additively. 

In summary, we have presented a simplified model for quo- 
rum sensing induced luminescence phenotypes in V. harveyi 
and V. cholerae. Our analysis suggests that a single parameter 
difference in our model effectively reproduces many features 
of observed luminescence curves which were hitherto uncon- 
nected. Thus large sequence-based differences are, in principle, 
not required to explain the dramatic differences between the lu- 
minescence phenotypes in these two species. Our model also 
makes testable predictions for observable luminescence pheno- 
types (specifically in V. cholerae) which, if validated, should 
shed new light on luminescence regulation by quorum sensing. 
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4. Appendix 

4.1. Single sRNA model 

Here we provide additional details for the single sRNA 
model. For convenience, we introduce the following dimen- 



= 1 - axy — X, 
= \-pxy-y. 



(7) 
(8) 



These equations can be readily solved to determine how steady 
state sRNA-mRNA levels change as system parameters are var- 
ied. In the limit ff,/? » 1, the solutions show a sharp transition 
as the ratio ajfi changes from (a/jS) < 1 to (ajfi) > 1. This 
parameter regime lets the system respond in an ultrasensitive 
manner as discussed in previous works I23I1 . During quo- 
rum sensing, the production rate of the sRNA (A: ^) decreases and 
hence the parameter /3 is lowered as bacteria make the transition 
from low-cell density to high-cell density. Correspondingly the 
system evolution traces out a trajectory in ia,/3) phase space. 
For Q',y6 » 1 the target mRNA levels show a sharp change as 
the line a = yS is crossed; thus it seems natural to identify the 
sharp transition observed in the luminescence profile with the 
sharp transition in target mRNA levels as /3 is lowered. How- 
ever, as argued in the previous sections, this identification is 
unlikely to be valid based on the following observations: 1) the 
quorum sensing response in V. harveyi is observed to be graded 
rather than all-or-none Issll . 2) Recent experiments have shown 
that HapR represses aphA at low cell density 111 ill , thus target 
mRNA levels are significant even at low cell densities. 3) north- 
ern blots show Uttle difference in the amount of sRNA in V. 
cholerae when the target mRNA (hapR) is deleted fTTj. 

Observation 3.) from above suggests that the sRNA and 
mRNA interactions occur in a parameter regime where the 
sRNA is never fully suppressed. To adhere to this constraint, 
we look at the limit a «; 1 and yS > 1 which effectively holds 
the sRNA concentration constant regardless of the target mRNA 
concentration. In this limit, the system no longer has an ultra- 
sensitive response, but instead responds in a controlled manner. 
As sRNA production rates double, the mRNA concentrations 
are about halved, thus allowing for a graded response. Apply- 
ing the limits a I and /? > 1 to equations Q and (jSj in the 
steady state explicitly shows the controlled response: 



X ^ 1, 

y ^ 1/(1 +yS). 



(9) 
(10) 



We note that recent research II28I1 has provided evidence for 
dosage compensation in V. cholerae due to regulation of sRNA 
production by HapR. While the inclusion of this effect will im- 
ply that the restrictions on the parameter a noted above are not 
required, the observation that the response of the regulated tar- 
get is graded in a controlled manner remains the same. 

4.2. Multiple sRNA with auto regulation model 

Here we provide additional details for the multiple sRNA 
with autoregulation model. In V. harveyi, there are a total of five 
sRNAs; however, only four are actively controlling the concen- 
tration of luxR mRNA. Likewise, V. cholerae contains four ac- 
tive sRNAs. Including multiple sRNAs has generated the new 
constants: k^., 7,, and yu.v,. However, we make the assumptions 
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that each sRNA has equal affinity to the target mRNA and all the 
sRNAs have the same degradation rate in both bacteria. These 
assumptions return y,- to y and n^, to jU j - To model autoreg- 
ulation, we introduce the dimensionless parameter yo as the 
threshold concentration for effective autoregulation of the tar- 
get gene. In dimensionless units, the system should be tuned 
in such a way that yb is not larger than the maximum value 
obtainable by y which is 1 . 

Using similar dimensionless parameters as the single sRNA 
model, we replace fi with - ('/k^.)/{fi^fiy), and introduce the 
dimensionless parameter e - fiy/fix, which is only necessary in 
the time dependent solutions of the model. Equations (|3]l and 
(lU are therefore rewritten as the following set of dimensionless 
equations 



LCD (1) 



HCD (2) 



Stationary (3) 



dxi 
z— = 1 
dt 

^ = _ 

dt 



■ axiy - Xi 
1 



1 + iy/yo) 



■y-y- 



(11) 

(12) 



The addition of multiple sRNAs to the model does not change 
the production rate of the target mRNA; therefore, we still con- 
sider the system to be in the parameter space where a « 1 . At 
steady state, i,- ^ 1 and the summation in equation (fTSl i reduces 
to Y^iPfy. Since y is independent of the summation, the sum 
is only of yS,, which results in just a constant representing all 
the contributions of the sRNAs, YiiPi ~^ Ptotai- The effects of 
multiple sRNAs are all integrated into the constant p total, and 
their removal via mutations to the wild-type bacteria is equiv- 
alent to reducing the maximum and minimum value of Ptotai as 
the bacteria moves from low cell density to high cell density 
respectively. The steady state concentration of y therefore be- 
comes: 



\ yoj l+Plotal 



(13) 



The effect of the autoregulation is best seen via different lim- 
iting cases of ratio ylyo in equation (fT3] ). When ^ 1, 
then ylyo » 1 resulting in y x ^jyoKl + Pwtai) ~ 0. This 
corresponds to the case where autoregulation is maximally on 
which prevents the system from sustaining any appreciable 
amount of protein. When yz) » 1, then ylyo *k 1 resulting in 
y K 1/(1 +j3,oiad which is similar in form to equation ( fTOl i where 
autoregulation is absent from the system. Since the amount of y 
is constrained to a value between and 1, andy/j is the effective 
concentration needed of the target protein before autoregulation 
occurs, we set yo - 0.75, which corresponds to the production 
rate dropping by close to half as seen by experiment js]] . 

4.3. Parameter space analysis 

Here we discuss the various parameter values and their as- 
sociated experimental motivation used in the preceding mod- 
els. Fluorescence experiments involving the expression of V. 
harveyi's qrr2 in the low-cell density and high-cell density lim- 
its provide a possible measure for estimating the sRNA fold 
change between the two cell density Umits ll33ll . The same type 
of fluorescence experiment shows the translational rate of luxR 
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Table 1: A table of the different /?,(/) values for WT, AluxU, AluxO, and the 
qrr mutants. 



113311 when LuxR autoregulation is removed. A direct determi- 
nation of relative fold differences using Real-Time Quantitative 
PCR for qrrl, qrr2, qrr3, qrr4, qrr5, and luxR with autoregu- 
lation intact has also been done jsoll . In the case without au- 
toregulation, luxR translational levels change ~10 fold and qrr2 
expression levels also change ~10 fold. 

With regards to the experimentally shown constraints, we let 
j6, change 10 fold between the low-cell density and high-cell 
density limits. Since we are in the limit where is always 
greater than 1, we chose y6, * 20 for the low-cell density result- 
ing in j6, a; 2 for the high-cell density limit. We set a = 0.1 to 
satisfy the previously discussed constraint: a <sc 1 . The values 
chosen for a and minimally satisfy the limits set on param- 
eter space; and yet, the system behaves in a manner consistent 
with experiment. Smaller values of a and/or larger values of jS, 
are also consistent with experiment showing robustness of the 
system in this parameter regime. 

To incorporate the effect of the system entering stationary 
phase, we introduced the parameter 6 such that k^. — > 6k ^. . we 
set 6 to the fixed value 0.025 - the maximal value necessary 
to have a clear enough distinction between the distributions at 
position (2) and (3) for the luxU mutant, where position (3) 
represents the colony entering stationary phase (see figure|6]l. 

a, Pi, and 6 are the only parameters necessary to determine 
the (normalized) mean values of LuxR/HapR. Furthermore, a 
and 5 remain a fixed value throughout our analysis, 0.1 and 
0.025 respectively. , which is a function of the sRNA produc- 
tion rates, only changes in value between the low-cell density 
limit and the high-cell density limit. The effects of the different 
mutants are also embedded into as they represent variations 
to the sRNA production rates relative to the WT. 

The critical factor in determining the decomposition of y6, is 
the fraction of LuxO (/) that is capable of promoting the pro- 
duction of sRNA. Therefore, y6, is a function of /, j8,(/). To bet- 
ter understand the contributions of LuxU and LuxO to y6, (/), we 
specify the different values /?,(/) can achieve depending on cell 
density and genotype. First, quantitative real-time PCR exper- 
iments show a basal rate sRNA production that is independent 
of the presence of LuxO which we label: jS,(0) jjoll . Next there 
is the rate, j8,(/o), that depends on the presence of LuxO which 
is evident in the AluxU mutant showing a wild-type like lumi- 
nescence phenotype in V. cholerae [21,] . Then there are the rates 
associated with phosphorylating LuxO, the dominant factor in 
sRNA production, in the low cell density limit (fiiifLco)) and in 
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the high cell density limit (fiiifHco))- The different values jS^/) 
for WT, AluxU, AluxO, and the qrr mutants are listed in table 1 . 
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